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", Abstract 

■ We present measurements of time-dependent asymmetries in decays of 

CsJ , neutral B mesons to the final states B° — > p^ir^. Measurements are based 

£^ | on a 140 /ft -1 data sample collected at the Y(4«S) resonance with the Belle 

detector at the KEKB asymmetric-energy e + e~ collider. 
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^ 1 Introduction 

in 

The angle 02 can be determined using a time dependent analysis in charmless decays 
such as B° — > 7r + 7r _ , B° — > p + 7r~ and 5° — »■ p + p~ ■ The extraction of 02 has been 
studied in several theoretical approaches. [TJ 12 El An isospin analysis that requires 
knowledge of the branching fractions of B° — > h°h°(h = it, p) can determine 02 with 
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relatively little theoretical uncertainty, 
i | Here we present a time dependent analysis of CP violation in B° — > p ± 7r =F decays. 

We also report on the first evidence for 5° — > p°7r° decay. 

: 2 £° -> p + tt- 

As the final state P° — > p + 7r~ is not a CP eigenstate, one has to consider four decay 
modes with different charge and flavor combinations: B° — > p^vr 111 and 5° — ► p^TrT. 
The decay rates can be written as 



7"*^ (At) = (l±Ag P 



e -|At|/r B o 



8r B o 

x{l - g(C p7r ± AC p7r ) coslAm^At) 
+q(S p7r ±AS pn )sm(Am d At)}, (1) 

where At = t p7T — t tag is the proper-time difference between the fully reconstructed 
and the associated B decay, and q = +1(— 1) corresponds to B°(B°) tags. The pa- 
rameters S p7T and C pw are associated with mixing- induced CP violation (related to 
2 ) and flavor- dependent direct CP violation, respectively. The parameters AS pn 
and AC pn are CP-conserving. AC p7T describes the asymmetry between the rates 
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T{B° -> P+7T-) + r(P° -> p-7r+) and r(P° -> p-7r+) + r(P° -> p+vr-). 
depends in addition on difference in strong phases between the amplitudes con- 
tributing to P — > pit decays. The direct CP violation is related to the time and 
flavor integrated charge asymmetry, Acp, defined as 

jV(p+7T-) - jV(p-7T + ) 

A7P ~ N(p+n-) + N(p-7;+y [ ] 

where N(p + 7t~) and N(p~7i + ) are the sum of the yields for P° and P° decays to p + 7r~ 
and p~7r + , respectively. We also define two separate direct-CP related asymmetries: 

iV(P° - p^) - iV(P° - p^) 



±_ ~ .Y(P" -> pT7r±) + /V(P° -> p : 77' 



1 ± ACp,,- i ^4-^-pCpT,- 



The analysis is based on the quasi-two-body approach. We consider B° — > 
(7r ± 7r°)7r T candidates from two distinct bands in the 7r + 7r~7r° Dalitz plot. The 
B° — > p ± 7r =F candidates are selected by requiring two oppositely charged pion- 
candidate tracks with 7r° candidates. Charged tracks are required to originate from 
the interaction point and have transverse momenta greater than 100 MeV/c. To 
identify tracks as charged pions, we combine information from time-of-flight coun- 
ters, aerogel threshold cherenkov counters and the central tracking chambers. Neu- 
tral pion candidates are reconstructed from photon pairs with invariant masses in 
the range 0.118 GeV/c 2 < M 77 < 0.150 GeV/c 2 , and momentum larger than 200 
MeV/c in the laboratory frame. Photon candidates are selected with a minimum 
energy requirement of 50 MeV in the barrel region of the electromagnetic calorimeter 
(32° < 6 1 < 129°) and 100 MeV in the endcap regions, where 1 denotes the polar 
angle of the photon with respect to the beam line. We also apply a cut on the an- 
gle between the photon flight direction and the boost direction from the laboratory 
system in the 7r° rest frame, | cos6^ ec | < 0.95. 

B meson candidates are identified using two kinematic variables: the beam 



energy constrained mass M bc = \J E£ eam — p B , and the energy difference AP = 
Eb — Pfeeam, where Eb eam is the center-of-mass (CM) beam energy, and Eb and ps 
are the reconstructed energy and momentum of the P candidate in the CM frame. 
We accept events in the region M bc > 5.21 GeV/c 2 and -0.3 GeV < AP < 0.2 GeV 
and define a signal region as M bc > 5.27GeV/c 2 and -0.1 GeV < AP < 0.08 GeV. 
To avoid the region where the p + ir~ and p~tt + amplitudes interfere, we exclude 
candidates with both M^± v o smaller than 1.22 GeV/c 2 . 

The dominant background to P° — > p + 7i~ is e + e~ — > qq (q = u, d, s, c) continuum 
events. To discriminate the signal from the background we use the event topology, 
which tends to be isotropic for PP events and jet-like for qq events. We combine a 
Fisher discriminant based on five modified Fox- Wolfram moments jl] with the cosine 
of the angle between the flight direction of the P in the CM frame and the electron 



2 



-0.3 -0.2 -0.1 0.1 -0.3 -0.2 -0.1 0.1 
AE (GeV) AE (GeV) 



0.2 




Figure 1: Projections for the results of the 2-D unbinned likelihood fit to the AE(top) 
and Mb c (bottom) distributions for p + n~ candidates (left) and p~n + candidates 
(right). 



beam direction to form the likelihood ratio 1Z = £ BB -/£ BB - + C q q. Here, C BB and 
Cgq are likelihood functions for signal and continuum events, respectively. The 
continuum background is reduced by applying a cut 1Z > 0.8. 

The 6-flavor of the accompanying B meson is identified using tracks not associ- 
ated with the reconstructed B° — > p ± ir T decays. The tagging algorithm is identical 
to that used for the Belle sin20i measurement jS] • The tagging information is rep- 
resented by two parameters, q — ±1, the flavor, and a quality factor r that ranges 
from r = for no flavor discriminant to r = 1 for unambiguous flavor assignment. 

The signal-side tracks and the tag-side tracks are fit for the signal (zcp) and 
tag (ztag) side vertices using an interaction point constraint. Taking into account 
that BB are produced practically at rest in the T(4S) CM frame, the proper time 
difference At is proportional to the distance between vertices Az = z cp — z tag '- At ~ 
Az/[3jc, where f3j = 0.425 is the Lorentz boost of the e + e~ system. 

The p ±r K T yields are extracted using an unbinned maximum-likelihood (ML) 
fit to the two-dimensional (Mb c , AE) distribution. Four types of background are 
considered: continuum qq, charm B background (b — > c), B — > pK, and rare, 
charmless B decays (rare B). The M^ c and AE distributions in the signal region 
along with projections of the unbinned ML fit are shown in Fig. ^ The fit yields 
483 ±46 B° — > p ± tt t candidates and a time and flavor integrated charge asymmetry 
A(J P = —0.16 ± O.lO(stat). Fig. |21 shows signal yields extracted from the fits to the 
(Mb c , AE) distributions for different M(7t ± 7t =f ) and cos6^ ei mass and helicity bins 
respectively. The result of B — > p ±r K T MC simulation are shown as histograms. 
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Figure 2: Signal yields as functions of M n ± n o (left) and cos6^ eZ (right). Histograms 
are expectations from the B — > p ± 7r =F MC simulation. 



The CP parameters C, AC, S, AS are obtained from an unbinned ML fit to 
the At distribution. The likelihood function is defined as 

N 
i=l 

+ f P KP p K{Atl) + f qq -Pgg(Ati) 

fsB^BB^^^i) frarePrarei^Ati)} , (4) 

where / ' pK qq - BB - rare is the fraction of events determined on an event-by-event 
DcLSlS cLS £L function of the M^ c and AE for each flavor tagging r interval and 
7T° momentum range. P pn (Ati) and P^ c (Atj) are the probability density func- 
tions (PDF) for signal reconstructed with right-sign p charge and wrong-sign p 
charge. They are determined as the convolution of the true PDF defined in Eq. 
with the At resolution function. 1 For P p ^(At) we assume C = S = AS = 0, 
AC = — 1, and A P q P = 0. The PDF for continuum background events is P qq -{At) = 
{(/r/2r 6A . g )e- |Ar|/Tbfc 9 + (1 - f T )S(At)}/2 convoluted with a resolution function R, 
where f T is the fraction of background with the effective lifetime Tbk g , and S is the 
Dirac delta function. The resolution function R for the background is parameter- 
ized as the sum of two Gaussian distributions. The parameters for the PDFs are 
determined from the sideband data for continuum and from the MC simulation for 
the BB and rare B decays. 

A fit to the 1215 candidates in the signal region (328.7 signal, 11.2 B — > pK, 
833 qq, 23.3 BB and 18.8 rare B events) yields C p7T = 0.25 ± 0.17^, A C p7T = 
0.38±0.18lg;^, S pn = -0.28 ±0.23^, and AS pn = -0.30 ± 0.24 ± 0.09. The first 
errors are statistical and the second systematic. The resulting At distributions are 
shown in Fig. El along with projections of the fit. 

1 the true PDF given by Eq.Q is modified to take into account the effect of incorrect flavor 
tagging 
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Figure 3: At distributions for B° — * p + ir~ (top) and B° — > p~ir + (second row), and 
resulting CP asymmetry (bottom rows). The CP asymmetry is shown separately 
for high and low quality flavor tags. The smooth curves are fit results. 



The systematic errors include the uncertainties in the background fraction, back- 
ground At PDF, wrong-tag fractions, vertex reconstruction, pir and pK At resolu- 
tion functions, and fitting bias. 

3 B°^p°7r 

The main reconstruction features of the analysis are the same as those for the 
B° —>■ p ± 7r =F decay channel. Neutral pion candidates are selected from the range 
0.115 GeV/c 2 < M 77 < 0.154 GeV/c 2 . Possible contributions from (b — > c) back- 
ground to the 7T + 7r _ 7r° final state are vetoed for the decays B° — ► D~tt + , D 7t° 
and J/tpir . We require the 7T + ir~ invariant mass to be in the range 0.5GeV/c 2 < 
771^+^- < 1.1 GeV/c 2 and the helicity angle for p° candidates | cos9^ el \ > 0.5, where 
Q p hel is defined as the angle between the 7r~ direction and the direction opposite the 
B in the p rest frame. The contributions from the B° — > p ± n T events are vetoed by 
rejecting candidates with 7r ± 7r invariant mass 0.5GeV/c 2 < m^i^o < 1.1 GeV/c 2 . 

To reduce the dominant continuum background we use the likelihood ratio TZ 
described in Sec. El Additional background discrimination is provided by the b- 
flavor tagging algorithm. We utilize the tagging parameter r, which ranges from 
to 1 and is a measure of the confidence that the b flavor of the accompanying B 
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Figure 4: Distribution of AE(M\ >C ) in the signal region of M^ C (AE). Projections of 
the unbinned ML fit are shown as the smooth curves; the dashed line represents the 
signal component; the dot-dashed line is the contribution from continuum events, 
and the dotted line is the composite of continuum and 5-related backgrounds. 



meson is correctly assigned. Events with a high value of r are well-tagged and are 
less likely to have originated from the continuum production. 

B candidates are accepted in the region 5.23GeV/c 2 < M bc < 5.3GeV/c 2 , 
-0.2 GeV < AE < 0.4 GeV. The signal region is defined as 5.269 GeV/c 2 < M bc < 
5.29GeV/c 2 and -0.135 GeV < AE < 0.082 GeV. The signal yield is obtained 
using an unbinned ML fit to the M^ C -AE distribution. The fit includes components 
for signal, continuum, b — > c decays and the charmless B decays B + — > p + p°, 
B + — > p + 7T° and B + — > 7r + 7r°. The fit result is presented in Fig. 0] The signal yield 
is 15 ± 4.8 events, corresponding to a 3.6a significance. 

We consider systematic uncertainties from the following sources: 7r° reconstruc- 
tion efficiency (±3.5%), track finding efficiency (±2.4%), not included charmless B 
decays (±5.3%), PDF shapes (±1.6%), and possible interference with B° — > p ± 7r T 
decays (9.3%). We obtain the branching fraction 

B(B° -> pV) = (5.1 ± 1.6{stat) 

±0.9(syst)) x 10~ 6 . (5) 



4 Summary 

We have measured CP violating parameters for B° — > p ± 7r T decays and have also 
obtained the first evidence for the decay B° — > p°7r°. The branching fraction mea- 
sured is relatively large (5.1±1.6±0.9) x 10 -6 , and thus the isospin analysis does not 
significantly limit the penguin diagram contribution. jHj A recent theoretical model 3J 
which assumes broken fiavor-5'C/(3) symmetry gives 02 — (102 ± 11 ± 15)°. The first 
error is experimental while the second is theoretical uncertainty. 



6 



References 

[1] J. Charles et al. hep-ph/04061841 N. G. Deshpande et al. Phys. Lett. B 473, 
141 (2000); A. Ali et al. Phys. Rev. D 59,(1999); M. Beneke and M. Neubert, 
Nucl. Phys. B 675, 333 (2003). 

[2] M. Gronau and D. London, Phys. Rev. Lett. 65, (1990); Y. Grossman and 
H. R. Quinn, Phys. Rev. D 58, (1998). 

[3] M. Gronau and J. Zupan, |hep-ph/040f002| 

[4] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 



[5] H. Kakuno et al, |hep^e x/0403022, 

[6] J. Stark, eConf C0304052, WG423 (2003). 



7 



